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Abstract 

We examine both the dynamical and the multifractal properties at the chaos thresh- 
old of logistic maps with general nonlinearity z > 1. First we determine analytically 
the sensitivity to initial conditions Then we consider a renormalization group 
(RG) operation on the partition function Z of the multifractal attractor that elimi- 
nates one half of the multifractal points each time it is applied. Invariance of Z fixes 
a length-scale transformation factor 2~ v in terms of the generalized dimensions Dp. 
There exists a gap Arj in the values of r\ equal to X q = 1/(1 — q) = — Dzlo 
where X q is the g-generalized Lyapunov exponent and q is the nonextensive entropic 
index. We provide an interpretation for this relationship - previously derived by 
Lyra and Tsallis - between dynamical and geometrical properties. 
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1 Introduction 



The thermodynamic framework for the characterization of multifractals [1] 
has proved to be a useful and insightful tool for the analysis of many complex 
system problems in which this geometrical feature plays a prominent role in 
establishing physical behavior. This has clearly been the case of dynamical 
systems, turbulence, growth models, sand pile models, etc. [2]. As it has re- 
cently been pointed out, multifractality is also of relevance to critical [3] and 
glassy dynamics [4]. In the former case the geometry of critical fluctuations is 
found to be multifractal and this has a hold on their dynamical properties [3], 
whereas in the latter case slow dynamics develops as phase-space accessibility 
is reduced and is finally confined to a multifractal subspace [4]. Concurrently, 
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in the study of the same types of systems the nonextensive generalization of 
the Boltzmann Gibbs (BG) statistical mechanics [5], [6] has recently raised 
much interest and provoked considerable debate [7] as to whether there is 
firm evidence for its applicability in circumstances where a system is out of 
the range of validity of the canonical BG theory. Slow dynamics is plausi- 
bly related to hindrance of movement in phase space that leads to ergodicity 
breakdown, and this in turn has been linked to the failure of BG statistics [4]. 
Here we briefly examine the relationship between multifractal properties and 
nonextensive dynamical evolution at the edge of chaos, a nonergodic state, in 
unimodal maps, as part of a wider analysis [8] , that might contribute to reveal 
the physical source of the nonextensive statistics. 



2 Nonextensivity of the Feigenbaum attractor 

An important class of multifractals are the strange attractors that are gener- 
ated by iteration of nonlinear maps. Amongst these, a prototypical example is 
the Feigenbaum attractor, also known as the edge of chaos, generated by the 
logistic map , / Mj2 (x) = 1 — fix 2 , — 1 < x < 1, when the control parameter \i 
takes the value ji^ = 1.40115... [1]. Recently, the predictions of the nonexten- 
sive theory have been rigorously proved [9] [10] for this state that marks the 
transition from periodic to chaotic motion and that is characterized by the 
vanishing of the ordinary Lyapunov exponent Ai = 0. At the edge of chaos 
(and also at infinitely many other critical states of the same map such as the 
pitchfork and tangent bifurcations where Ai = [11]) the sensitivity to initial 
conditions £ t = \dx t /dx \ (where x t is the orbit position at time t given the 
initial position xq at time t — 0) generally acquires power law instead of ex- 
ponential time dependence. The nonextensive formalism suggests [12] that ^ 
when Ai = obeys the g-exponential expression, 

6 = exp q (X q t) = [1 - (q - l)A^]- 1/(9-1) , (1) 

containing the entropic index q and the g-generalized Lyapunov coefficient X q . 
The BG exponential form £ t = exp(Aii) is recovered when q — > 1. Further, it 
was conjectured [12] that the ordinary Pesin identity X± — K ± , where Ki is the 
Kolmogorov-Sinai (KS) rate of entropy increase [13] is replaced by \ q = K q , 
where K q t = S q (t) — S q (0), t large, and where 

^E«'n,(i)=i^M (2) 

is the Tsallis entropy (ln g y = {y l ~ q — 1)/(1 — q) is the inverse of exp (?/)). 
In the limit q -> 1 K g becomes K x = f -1 [Si(t) - Si(0)] where S^t) = 



2 



— Y^Li Pi(t) lnpi(t). In Eq. (2) p^t) is the probability distribution obtained 
from the relative frequencies with which the positions of an ensemble of trajec- 
tories occur within cells % = 1, W at iteration time t. Based on the renormal- 
ization group (RG) transformation R/(x) = af(f(x/a)) where a ~ 2.5029 
is the Feigenbaum's universal constant that measures the power-law period- 
doubling spreading of iterate positions, a rigorous analytical confirmation of 
Eq. (1) and also of X q = K q have been obtained [9] [10]. Specifically, \ q and q 
are simply given by X q = In a/ In 2 ~ 1.3236 and q = 1 — In 2/ In a ~ 0.2445. 

On the other hand, a remarkable link between the well-known [1] geometrical 
description of multifractals and the nonextensive statistics has been derived 
and corroborated numerically not only for the Feigenbaum attractor but for 
the edge of chaos of the generalization of the logistic map to nonlinearity 
z > 1, fn, z (x) — 1 — fi\x\ z , —1 < x < 1 [14]. This consists of the expression 
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where and -D-oo are the extreme values of the generalized dimension 
spectrum Dp, and correspond to the smallest and largest length scales of 
the multifractal set [1]. For the z- logistic map one has = ln2/zlno; and 
-D_oo = ln2/lno; (notice that a depends on z) [1]. In view of the clear asso- 
ciation of the nonextensive statistics with the properties at the edge of chaos 
of the logistic map we have pursued the study of this critical state by means 
of the standard thermodynamic tools that characterize multifractals. Here we 
advance some of the results we have found [81. 



3 Sensitivity to initial conditions 



We determine the index q and the g-Lyapunov coefficient X q at the edge of 
chaos of the z-logistic maps through direct evaluation of the sensitivity to ini- 
tial conditions £ t . By making use of the properties [15] g(0) = (— a) n g^ 2n \Q) = 
1, g{\) = (- a ) n g {2n \a- n ) = -a~ l and 

g\l) = {-lYg'{a- n )g'{g{a- n ))...g'{g^- 1 \a- n )) = -a*~\ (4) 



of the fixed-point map, g(x) = lim n ^ 00 (— a) n fjf™\— x/a n ), into the definition 

of 6, 



6 = 



dx t 




dg^-^ix) 


dx 




dx 



(5) 

x=xq 
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Fig. 1. a) Absolute values of \xt\ vs t in logarithmic scales for the orbit with initial 
condition xo = at fj,^ of the logistic map z = 2. The labels indicate iteration 
time t. b) Same as a) with \xt\ replaced by 1 — \xt\- 

with the choice x = 0, and for iteration times of the form t — 2 n — 1, we 
obtain £ t = a^ z ~ V)n . This result can be re-expressed as 



& = exp,(v) = [i-(? -i)vr 



-V(9-l) 



(6) 



with 



1-9 



= (z-l) 



In a 
TnT' 



(7) 



For details see Ref. [8]. Thus, without reference to the multifractal properties 
of the attractor, we have corroborated the result in Eq. (3) [14] and obtained 
additionally the g-Lyapunov coefficient X q — (1 — Alternatively, £ t can 

be obtained from |Ax t | / \Axq\ in the limit A^o — > where Axj is the distance 
between two orbits at time t when they were a distance Axo ait — 0. Following 
the procedure described in Ref. [9] with the choices \Axq\ = a - - 7 ' — a~ l , i,j>0 
and t = 2 n - 1, n = 0, 1, we obtain \Ax t \ = \Ax \ a {z ~ 1)n [8]. Next we study 
a multifractal property of the attractor and make contact again with Eq. (3). 



4 Scaling properties at the edge of chaos 



In Fig. la we have plotted (in logarithmic scales) the first few (absolute values) 
of iterated positions \xt \ of the orbit at //qo starting at xq = where the labels 
indicate iteration time t. Notice the structure of horizontal bands, and that, in 
the top band lie half of the attractor positions (odd times), in the second band 
a quarter of the attractor positions, and so on. The top band is eliminated by 
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functional composition of the original map, that is by considering the orbit 
generated by the map g^ 2 \0) instead of <?(0). Successive bands are eliminated 
by considering the orbits of g^ 2 \0), k = 1,2.... The top band positions (odd 
times) can be reproduced approximately by the positions of the band below 
it (times of the form t — 2 + 4n, n — 0, 1,2...) by multiplication of a factor 
equal to a, e.g. \x\\ ~ a ja^l- Likewise, the positions of the second band are 
reproduced by the positions of the third band under multiplication by a, e.g. 
| #2 1 — ot |^4 1, and so on. In Fig. lb we show a logarithmic scale plot of 1— 
\x t \ that displays a band structure similar to that in Fig. la, in the top band 
lie again half of the attractor positions (even times) and below the other half 
(odd times) is distributed in the subsequent bands. This time the positions in 
one band are reproduced approximately from the positions of the band lying 
below it by multiplication of a factor a z , e.g. 1 — \x 3 \ ~ a z (l—\x 5 \). The scaling 
properties amongst iterate positions merely follow from repeated composition 
of the map (for more details see Ref. [8]) and suggest the construction of 
an RG transformation described below. For this purpose it is convenient to 
write the two scale factors a and a z as 2 v -°° and 2 Vcc , respectively, where 
i]_ 00 = In a/ In 2 and r]^ = zlna/ In 2. 



5 'Take away half RG transformation 

Consider the positions of the orbit at //qo starting at xo = up to time 
t = 2 N — 1 and from them extract the set of lengths lf\ j = l...,2 N , that 
join closest adjacent positions. Each position is used only once so that the set 
covers the attractor but not all the interval — 1 < x < 1. To each length 

if^ assign the equiprobability = 2~ N . Following the standard procedure 
[1] we define the partition function 

4 0) (Ar) = £[pfn/fr, (8) 

and determine the (finite size) generalized dimensions Dp^ = t/(/3 — 1) by 
requiring Z^\(3,r) = 1, {Dp = limjy^oo Dp N ). In Fig. 2 we show results for 
D/3 t N for several values of z. 

We introduce now repeated length scale and restoring probability transforma- 
tions 

if =2~ k nf and pf) = 2~ k pf\ k=l...,2 N , (9) 

on the partition function Z^\f3,r). That is, to transform Z^\(3,t) into 
Z$({3, t) replace if and pf in Eq. (8) by lf ] = 2~nf ) and pf = 2- x pf , re- 
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Fig. 2. The generalized dimensions Dp jN for z = 1.25 (a(1.25) = 4.8323, N = 13), 
z = 1.5 (a(1.5) = 3.4479, N = 16) and z = 2 (q(2) = 2.5029, N = 14). 
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Fig. 3. The length rescaling index 77 in Eq. (10) for z 
(a(1.5) = 3.4479, N = 16). 



1.5 



spectively, and proceed similarly any number k of times. In the limit N — > oo 
the transformation eliminates one half of the multifractal points each time it 
is applied. By requiring Z^\(3,t) = 1 we obtain, again when N — ■> oo, the 
following result for the rescaling exponent rj of the interval lengths 



T 



(10) 



As shown in Fig. 3, i] approaches two limiting values, 77-00 = Dz\o and 77^ = 
as (3 — > — oo and (3 — > oo, respectively, and exhibits a singularity at /3 — 1. 
Except for values of /3 in the vicinity of /3 = 1 the exponent ?^ is always close 
to either i]-^ or 7]^ which correspond to the scaling factors lj ^ = cT k lf^ and 

l^P = a~ zk lf^ described in the previous section. As also shown in Fig. 3, there 
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Fig. 4. Successively RG transformed dimensions for k/N = 1 and z = 1.5. 



appears a gap in the permissible values of rj, 

1 1 

A77 = Tfo - V-oo = 



D. 



that is noticeably identical to the dynamical result Eq. (7), i.e. \ q = (1—q) 1 = 
Ar). 



Repeated application of the RG transformation when the limit N - 
taken before k — > oo leaves the generalized dimensions invariant, i.e. Dp — ^ p 
for all k. However when both k, N — > oo with k = N the RG transformation 
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leads to a fixed point — > Dp"' that has the form of a step function, 

D { ™ ] = In 2/2 In a for (3 < 1 and D { ™ ] = \n2/2z\na for (3 > 1, with a 
divergence at f3 — 1. See Fig. 4. [8]. 



6 Summary and discussion 



Assisted by the known properties of the fixed-point map g(x) we derived the 
expression for the sensitivity to initial conditions for the z-logistic map at 
the edge of chaos. As a consequence we corroborated the known expression 
for the entropic index q in terms of a(z) [14] and obtained in addition the 
closely related expression for the g-Lyapunov coefficient \ q (known previously 
only for the particular case z — 2 [9] [10]). Our results are rigorous. The 
expression obtained is precisely that first derived, heuristically, by Lyra and 
Tsallis [14], Eq. (3), that relates q to the smallest and largest -D_oo gen- 
eralized dimensions of a multifractal set. In order to probe this connection 
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between dynamic and geometric properties of the multifractal attractor we 
introduced an RG approach that consists of sequentially removing remaining 
halves 2~ k , k = 1,2, ... of points of the multifractal set while keeping the value 
of its partition function Z^\(3,t) fixed. We determined the length rescaling 
index 7] as a function of (3 and observed that there is a gap Arj of permissi- 
ble values in this quantity equal to X q — 1/(1 — q), as the edges of the gap 
correspond to r^-oo = Dz^ and r]^ = D^. This result can perhaps be better 
appreciated if we note that the orbit positions \x t \ evolve according to two 
basic alternating movements observable at time intervals of the form 2 n , one 
drives the iterate towards \x\ = 1 as a~ zn = 2~ Vac n and the other one towards 
x = as a~ n = 2~ r? -°° n . See the main diagonal subsequence positions in Figs, 
la and lb. The net growth in the distance between two orbits is the ratio of 
these two power laws, that is |Ax t | = |Ax | 2 Av n (see Section 3). Thus, Eq. 
(3) is obtained as the result of orbits that successively separate and converge, 
respectively, in relation to D_ OQ and D^. The power law for £ f can only persist 

if D_oo ^ D 

OO • 
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